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ORIGINAL ARTICLE 

ASXLl and DNMT3A mutation in a cytogenetically 
normal B3 thymoma 

R Belani\ G Oliveira^ GA Erikson^ S Ra^ MS Schechter^ JK Lee^ WJ Shipman^ SM Haaser^ and ATorkamanP 

The molecular drivers of thymoma are poorly understood. Outside of the identification of rarely occurring epidermal growth factor 
receptor and v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog mutations via candidate gene sequencing, 
mutations in common cancer genes have yet to be observed. Only a single thymoma genome sequence has been previously 
reported, with no mutations in known cancer genes identified. Thus, we attempted to identify somatic driver mutations in a 
cytogenetically normal thymoma. A stage IVB type B3 thymoma from a 47-year-old male of Asian descent with no history of 
myasthenia gravis or other autoimmune condition was genomically evaluated. Exome sequencing and low-pass whole-genome 
sequencing was performed to identify somatic point mutations, copy number changes and structural variants. Mutations in known 
tumor suppressors DNMT3A (p.G728D) and ASXLl (p.E657fs), consistent with mutations of known consequence in acute myeloid 
leukemia, were identified. Contrary to a previous report, this finding suggests the genetic etiology of thymomas may not be 
fundamentally distinct from other tumor types. Rather, these findings suggest that further sequencing of cytogenetically normal 
thymoma samples should reveal the specific molecular drivers of thymoma. 

Oncogenesis (2014) 3, elll; doi:1 0.1 038/oncsis.201 4.25; published online 7 July 2014 



INTRODUCTION 

Thymomas are rare tumors arising from the epithelial cells of the 
thymus, commonly associated with myasthenia gravis or other 
autoimmune disorders.^ The WHO classification system 
categorizes thymomas (type A, 81, 82, 83 and A8) and thymic 
carcinoma (type C) according to morphology, degree of atypia, 
histological similarity to normal thymic cells and degree of 
lymphocyte infiltration.^ Thymomas are staged (stage I, IIA, 118, 
III, IVA and IV8) according to degree of encapsulation and 
anatomic extent of disease with encapsulation decreasing and 
dispersion from the thymus increasing with stage.^ Late stage 
thymomas are very rare and tend to locally recur rather than 
metastasize. Type 83 thymomas are the most aggressive of the 
thymomas. 

Surgical resection is the favored treatment for thymoma. 
Complete surgical resection can be achieved for the vast majority 
of stage I and II thymomas. Stage lll-IV thymomas that cannot be 
fully resected are generally treated with radiation and platinum- 
based combination chemotherapy regimens leading to substantial 
tumor shrinkage;^"^ however, similar survival is observed with pre- 
operative radiation therapy without chemotherapy.^'^ Currently, 
there are no good biomarkers for response to therapy. Amplifica- 
tions and rare activating mutations of EGFR and overexpression of 
HER2, KIT, BCL2 and TP53, as well as loss of CDKN2A have been 
previously described.^'^° However, no strong evidence supports 
response of thymomas to EGFR or KIT targeted therapies." 

Thymomas and thymic carcinomas have been extensively 
profiled for chromosomal abnormalities, especially via 
comparative genomic hybridization,^^"^"^ revealing that the vast 
majority of 83 thymomas display chromosomal abnormalities. 
However, there is no information on tumorigenic point mutations 



in 83 thymomas.^ Only a single genome sequence of a 83 
thymoma has been reported thus far, for a tumor that displayed 
the typical chromosomal abnormalities for this tumor type (that is, 
gain of chromosome 1q, loss of chromosome 6 and 13q, and so 
on).^^ No mutations in common cancer-associated genes were 
observed, leading the authors to suggest that thymomas may 
evolve through distinct molecular mechanisms as compared with 
other tumor types. In this article, we report the identification of 
DNMT3A and ASXLl mutations in a cytogenetically normal stage 
IV8 type 83 thymoma. This is the first identification of driver point 
mutations in 83 thymoma and suggests that further genomic 
profiling of cytogenetically normal 83 thymomas may reveal the 
specific molecular drivers of this tumor type. 



RESULTS 

The patient presented with a 3-month history of cough. He did not 
have shortness of breath, hemoptysis or weight loss at presenta- 
tion. He also did not have neurologic symptoms consistent with 
myasthenia gravis. Family history is significant for thymoma in the 
mother. A computed tomography of the chest showed a 
mediastinal mass measuring 10.5 x 7.8 x 12.1 cm and a large 
(12 X 9.5 X 8.3 cm) left renal mass (Figures la and lb). The patient 
underwent computed tomography-guided biopsy of the mediast- 
inal mass that revealed mildly atypical, uniform spindled cells 
associated with dense collagenous tissue. The tumor cells showed 
immunohistochemical positivity for AE1/AE3 and cytokeratin-7. 
Immunohistochemical markers for renal cell carcinoma were 
negative. The pathologic findings were consistent with thymoma. 
Subsequently, the patient underwent removal of mediastinal mass 
and right upper lobectomy as an en bloc resection. There was a 
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separate right diaphragmatic nodule that was also excised. At the 
time of surgery, it was noted that the tumor was densely adherent 
to the pericardium and invaded the upper lobe of the right lung. 
Microscopic examination showed spindled to polygonal cells in 
varying sized lobules with accompanying fibrous bands of stroma 
invading the pericardium, pleura and lung (Figures 1c and d). Only 
occasional interspersed lymphocytes were seen. Focal lymph- 
vascular invasion was identified. The tumor cells showed ovoid 
nuclei with mild pleomorphism and scattered mitoses. The final 
margins were clear. The diaphragmatic nodule demonstrated 
tumor cells morphologically similar to the primary tumor. The final 
Masaoka stage of IVB was based on the finding of a separate 
distinct diaphragmatic tumor nodule and pleural invasion. The 
tumor had a WHO histologic category of B3 based on the 
predominance of epithelial cells and paucity of intraepithelial 
lymphocytes. No treatment was given before genome sequencing. 

Tumor genome overview 

An overview of the somatic variants observed is presented in 
Figure 1e. Overall, 14 somatic indels and 176 somatic single- 
nucleotide polymorphisms were identified, 17 of which impacted 
the amino-acid sequence of proteins (Table 1). Sixteen copy 
number gains, 1 copy number loss (Table 2) and 6 structural 
variants (Table 3) were also observed. Sequencing of the renal 
mass revealed no overlap in somatic mutations, confirming the 
tumors were of independent origin. No major interchromosomal 
events were present in the thymoma genome, nor were any major 
chromosomal aneuploidies frequently observed in B3 thymo- 
mas.^^ Structural variants that were detected did not contain 
breakpoints within genes of known relevance to tumorigenesis 
(Table 3). Thus, somatic events of interest were comprised of point 
mutations and subchromosomal copy number variants. 



Coding mutations 

Fourteen nonsynonymous, 1 nonsense and 2 frameshift somatic 
mutations were observed in the tumor (Table 1). Of these 
mutations, four stood out with clear functional relevance — 
frameshift mutations in FAI\/n93A and ASXLl, a nonsense mutation 
in GRMS, and a nonsynonymous mutation in DNMT3A. Although 
the mutations in FAIV1193A and CRMS clearly impact protein 
function, the relevance of these mutations to tumorigenesis is not 
immediately apparent. On the other hand, the somatic mutations 
in DNMTSA and ASXLl are of the type previously observed in other 
tumor types and expected to contribute to tumorigenesis. 

DNI\/1T3A is a DNA methyltransferase required for de novo 
methylation and mammalian development.^^ It is known to act as 
a tumor suppressor in acute myeloid leukemia (AML) and is 
observed mutated less frequently in other tumor types.^^'^^ The 
nonsynonymous mutation observed in this thymoma sample 
leads to a nonconservative substitution of aspartate for glycine at 
amino-acid 728 (p.G728D) within the DNA methylase domain of 
DNIV1T3A. This region is known to be responsible for the 
hydrophobic interaction between DNIV1T3A and its regulatory 
subunit DNIV1T3L, an interaction that is required for activation of 
DNIV1T3A methyltransferase activity and previously shown to be 
disrupted by mutations in this region.^ ^ Moreover, previously 
observed nonsynonymous mutations in AML cluster tightly 
within this region (shown in Figure 2)^°"^^. Strikingly and 
unusually, this somatic mutation was observed in homozygous 
state in the tumor: 61 of 62 (98.4%) exome sequencing reads 
from the tumor supported the mutation, 0 of 41 (0%) of 
exome sequencing reads from the normal sample supported 
the mutation, the homozygous state was confirmed via 
Sanger sequencing (Figure 2, inset), and no evidence for loss 
of heterozygosity was observed in the exome or low-pass 



Table 1. Coding mutations 



Gene 


Gene impact 


Protein location 


Prediction of impact 


Coordinates (hg19) 


ALDH1A3 


Missense 


T51I 


SIFT Condel 


chr15: 101 425 524 snp C>T 


AP3M2 


Missense 


H19N 


Polyphen, SIFT, Condel 


chr8: 42 012 260 snp C>A 


ASXLl 


Frameshift 


E657fs 


Damaging 


chr20: 31 022486 del G>- 


ATRIP 


Missense 


G21R 


SIFT Condel 


chr3: 48488310 snp G>C 


DNAJB12 


Missense 


L292V 


Neutral 


chrlO: 74098021 snp G>C 


DNMT3A 


Missense 


G728D 


Polyphen, SIFT, Condel 


chr2: 25 463 310 snp C>T 


FAM193A 


Frameshift 


H303fs 


Damaging 


chr4: 2 648428 ins ->T 


GRMS 


Nonsense 


R351* 


Damaging 


chrll: 88 386432 snp G>A 


HECTD3 


Missense 


R14P 


SIFT Condel 


chrl: 45 476889 snp C>G 


JPH2 


Missense 


V239M 


Neutral 


chr20: 42 788 712 snp C>T 


MAGI2 


Missense 


R375Q 


Condel 


chr7: 77 807 386 snp C>T 


M0RC1 


Missense 


R597H 


Neutral 


chr3: 108 724077 snp C>T 


NCALD 


Missense 


R181C 


Polyphen, SIFT, Condel 


chr8: 102 701 578 snp G>A 


PVRL1 


Missense 


L360V 


Neutral 


chrll: 119510648 snp G>C 


SH3BP1 


Missense 


T626P 


Neutral 


chr22: 38051 466 snp A>C 


TBC1D8 


Missense 


R144G 


Neutral 


chr2: 101 670 771 snpT>C 


TBKBP1 


Missense 


K75T 


Polyphen, SIFT, Condel 


chrl 7: 45 773 702 snp A>C 



Somatic protein coding mutations. 



< 

Figure 1. (a) Thymoma imaging, fluorodeoxyglucose positron emission tomography/computed tomography (CT) demonstrates intensely 
hypermetabolic anterior mediastinal mass on fused coronal and axial images (arrowheads). Cystic tumor component on attenuation- 
correction CT image is non-hypermetabolic on FDG PET (arrows), (b) Renal tumor imaging. CT scan with contrast (left) demonstrates large 
contrast-enhancing left renal mass (arrowheads) with posterior hydronephrotic upper pole collecting system (block arrow). Fused axial FDG 
PET/CT (right) demonstrates mildly hypermetabolic posterior tumor component (arrow). Upper pole hydronephrosis contains excreted urinary 
FDG (block arrow), (c, d) Thymoma histology, (c) Lobules of thymoma tumor cells infiltrate into the adjacent lung parenchyma (arrows: 
hematoxylin and eosin (H&E): x 5 magnification), (d) Sheets of thymoma type 83 polygonal tumor cells with ovoid nuclei (H&E: x 20 
magnification), (e) Tumor genome overview. Circos plot overview of somatic mutations identified in the thymoma. Outer text denotes 
location and gene impacted by somatic protein coding mutations. Middle gray ring depicts copy number variation (CNV) status (black = 2 
copies, red = CNV gains, blue = CNV losses). Inner text displays genes overlapping identified CNVs. Inner red lines display identified structural 
variants. 
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Table 2. Copy number variants 



Type 


Coordinates 


Size (Kb) 


Dosage 


Overlapping genes 


Gain 


cnri: 4/uuui-duuuuu 


1 30 


5 


None 


Gain 


cn rz. 1 4d 1 oU UU l - l 4d o l U UUU 


1 on 
1 oU 


o 
D 


Z.tDZ 


Gain 


cnrz. 1 DO zdu UU 1 - 1 Do 3/u UUU 


1 lU 


•3 
D 


LY 1 Ir 


Gain 


cnrz. 1 yz / /u UU 1 — 1 yz y4u UUU 


1 /U 


■3 
D 


1 IVItrrz 


Gain 


cnrz. z4U 1 1 U UU 1 -z4U z4U UUU 


1 dU 


•3 
D 


nUAL4 


Gain 


cn r4. 1 DO zzU UU l - l do ozU UUU 


1 UU 


3 
D 


rDAVV/ 


Gain 


cnr/: z/ izuuu i-z/ z4uuuu 


120 


3 


huaA gene cluster, iviiKiyDD 


Gain 


cnr/. DO zyu UU 1 -DO oyu UUU 


1 UU 


3 
D 


TA DD 

1 AKr 


Gain 


chr8: 142 200001-142 320000 


120 


3 


SLC45A4, DENND3 


Gain 


chr9: 75 690001-75 830000 


140 


3 


ANXA1, ALDH1A1 


Gain 


chr9: 139420001-139520000 


100 


3 


N0TCH1, MIR4674 


Gain 


chr13: 114 740 001-114920000 


180 


3 


RASA3 


Gain 


chr15: 99410001-99510000 


100 


3 


IGF1R 


Gain 


chr17: 80 740001-80920000 


180 


3 


TBCD, ZNF750, B3GNTL1 


Loss 


chr19: 54070001-54200000 


130 


1 


ZNF331, DPRX, I\/1IR520 cluster 


Gain 


chr21: 44690001-44900000 


210 


3 


SIK1 


Gain 


chr22: 46380001-46550000 


170 


3 


MIRLET7B, MIR3619, MIR4763 



Tumor-specific copy number variants, dosage and contained genes. 



Table 3. Structural variants 


Type 


State 


Breal<point 1 


Breal<point 2 


Overlapping genes 


Coamplicon 


Bal 


chr16: 26382 735-26386932 


chr16: 26 382 736-26386932 


None 


Inversion 


Unbal 


chr2: 87 651 839-87 656090 


chr2: 90482 288-90485 529 


Many 


Inversion 


Unbal 


chr2: 33141 477-33141 679 


chr2: 242 81 9 853-242 820 060 


Many 


Coamplicon 


Bal 


chr7: 158 706424-158 710034 


chr7: 158 706425-158 710035 


WDR60 


Insertion 


Bal 


chr7: 158 704 999-158 708434 


chr7: 158 705 000-158 709994 


WDR60 


Insertion 


Bal 


chr9: 139994 967-139997 831 


chr9: 139 994968-139 998202 


MAN1B1 


Abbreviations: Bal, balanced; Unbal, unbalanced. Tumor-specific structural variants. 



DNMT3A 



chr2:25463310C>T 



Tumor 



G728D 



vMvWliWkMM; 


Normal ' 

• m ]] 




niTT! 


I 




RAAA/P 




DNA_methylase 





400 
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800 



912 aa 



Figure 2. The DNMT3A p.G728D nnutation observed in this thynnoma sample is visualized in the context of other DNMT3A nnutations observed 
in AML genome sequences from TCGA. The protein sequence and functional domains are depicted on the x axis. The number of AML 
mutations is depicted on the y axis. Red circles correspond to truncating nnutations. Green circles correspond to nnissense mutations. Circle 
height corresponds to the number of mutations per position, however, the G728D indicator (black) is only nneant to indicate position of this 
mutation. Note the clustering of AML nonsynonynnous mutations around position 728. Mutations in this region reduce DNMT3A activity by 
disrupting the interaction between DNMT3L and DNMT3A. The Sanger sequencing validation trace of p.G7238D is also shown, demonstrating 
validation of p.G7238D as a homozygous somatic mutation. 



whole-genome sequencing data. Moreover, no other somatic 
mutations were observed at >60% allelic frequency in the tumor, 
suggesting a specific and early gene conversion event at this 
locus. Finally, a previously sequenced AML, reported by Hou 
et al.,^^ carried a mutation at this position (p.G728R) and at a 
nearby position (p.F731 L), supporting the notion that homozygous 
or compound heterozygous mutations at the DNMT3A — DNMT3L 
interface are required to effectively disrupt methyltransferase 
activity.^^ Thus, the type, location and unusual zygosity of this 
somatic mutation strongly suggest mutations in DNMTSA have a 
causal role in the etiology of thymoma. 



ASXL1 is a chromatin binding Polycomb group protein involved 
in transcriptional repression of numerous cell-type-specific sys- 
tems, most likely through the recruitment of methylated histone 
H3 to promoters of target genes.^^ ASXLl is a known tumor 
suppressor in myelodysplastic syndrome and chronic myelo- 
monocytic leukemia although it is frequently mutated in other 
tumor types such as lung squamous and bladder cancer.^^'^^ The 
heterozygous deletion of a single guanine nucleotide within exon 
12 of ASXLl leads to out-of-frame translation beginning at codon 
657 and premature truncation of the protein, removing >50% of 
the protein including the PHD domain responsible from histone 
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chr20:31022486 G>- 
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Figure 3. The ASXL1 p.E657fs nnutation observed in this thynnonna sannple is visualized in the context of other ASXL1 nnutations observed in all 
tunnor genonne sequences catalogued in TCGA. The protein sequence and functional donnains are depicted on the x axis. The nunnber of TCGA 
nnutations is depicted on the y axis. Red circles correspond to truncating nnutations. Green circles correspond to nonsynonynnous nnutations. 
Purple circles correspond to nnutations that are both nonsynonynnous and truncating in different gene isofornns. Circle height corresponds to 
the nunnber of nnutations per position, however, the E657fs indicator (black) is only nneant to indicate position of this nnutation. Note the 
clustering of TCGA truncating nnutations around position 657. The Sanger sequencing validation trace of p.E657fs is also shown, 
dennonstrating validation of p.E657fs as a heterozygous sonnatic nnutation. 



interactions. This mutation was confirmed via Sanger sequencing 
(Figure 3, inset). Similar heterozygous C-terminal truncations 
are sufficient to induce myelodysplastic syndrome.^^ Previously 
observed frameshift or nonsense mutations in ASXL1 across 
numerous tumor types catalogued in The Cancer Genome Atlas 
(TCGA) have been observed both upstream and downstream of 
amino-acid 657, with the highest concentration centered at 
position 657 (Figure 3). In fact, the majority of previously observed 
ASXLl mutations occur in the same exon as this mutation, exon 
^2 26,27 -p|^^3^ |3Q|.^ ^ypg location of this somatic mutation 
strongly suggests mutations in ASXLl have a causal role in the 
etiology of thymomas. 

Copy number variants 

A number of focal copy number gains were observed in this 
tumor, the most interesting being an amplification of the HOXA 
cluster (see Discussion section).^^ Although a number of genes are 
known to be involved in tumorigenic processes are observed in 
these focal amplifications, no significance can be confidently 
assigned to these events. 



DISCUSSION 

Our results show that mutations in DNMT3A and ASXLl have a role 
in the development of thymoma. Our findings contrast starkly 
with a previous genome sequence from a B3 thymoma.^ ^ 
However, the previously reported thymoma displayed 
cytogenetic abnormalities commonly observed in B3 thymomas, 
whereas this tumor displayed largely normal cytogenetics. These 
results are reminiscent of the preferential occurrence of driver 
point mutations in cytogenetically normal AML and suggests 
further sequencing of cytogenetically normal thymomas are fertile 
ground for the identification of the molecular drivers of thymic 
cancer.^^ No obvious germline variant predisposing the patient to 
thymoma or renal cell cancer was observed, nor was there any 
overlap in somatic mutations between the thymoma and renal cell 
cancer to suggest shared etiology of these two tumors. 

The first confidently identified driver mutation was a homo- 
zygous nonsynonymous mutation in DNMT3A. DNIV1T3A is a DNA 
methyltransferase essential for development, especially for 
de novo DNA methylation in stem cells.^^ DNMT3A is mutated in 
numerous tumor types, but is mutated especially frequently 
in AIV1L.^°'^^ DNIV1T3A and its binding partner, DNMT3L form 
a tetramer: DNIV1T3L:DNIV1T3A:DNIV1T3A:DNMT3L.^^ The most 
common DNIV1T3A mutations in AML occur at p.R882. Mutations 
at this position disrupt binding at the DNMT3A:DNMT3A interface, 
and impact activity in dominant-negative manner,^^'^° or at least 
reduce the odds of proper tetramer assembly by ~75% as a 



single mutated DNMT3A protein is sufficient to disrupt tetramer 
assembly. Alternatively, the DNMT3L:DNMT3A interface is 
disrupted by somatic mutations within the mutational hotspot 
region centered on the p.G728D mutation observed in this 
thymoma.^^ Heterozygous mutations at this site may not be 
expected to influence activity as markedly as mutations at p.R882 
as they likely allow for the formation of partially complete 
complexes, for example, DNMT3A:DNMT3A:DNMT3L complexes 
could form 50% of the time. Thus, it is notable that the p.G728D 
mutation was observed in homozygous state in this tumor. 
Homozygous somatic mutations are highly unusual, suggesting 
DNA repair via double-stranded break repair or synthesis- 
dependent strand annealing and/or some other gene 
conversion event occurred specifically at this site given we did 
not find any evidence for loss of heterozygosity. 

The second confidently identified driver mutation was a 
frameshift mutation in exon 12 of ASXLl. ASXL1 is a known tumor 
suppressor in hematological malignancies, where the majority of 
frameshift and nonsense mutations are also observed in exon 
^2 26,27 ^5XL1 loss leads to loss of H3K27me3 and increased 
expression of genes repressed by these histone marks, especially 
the anterior HOXA gene cluster.^^ HOXA genes themselves are 
capable of leukemogenic transformation.^^ Interestingly, we also 
observe a concomitant and specific amplification of the HOXA 
gene cluster in this tumor. 

It is intriguing that both DNMT3A and ASXLl are known tumor 
suppressors in hematological malignancies^^ and that the 
observed mutation characteristics are consistent with those seen 
in AML. However, both genes have been observed mutated, albeit 
at lower frequencies, in epithelial malignancies as well. Although 
co-occurrence of DNMJ3A and ASXLl mutations are quite rare in 
epithelial and hematological malignancies, sequencing of a larger 
number of thymomas will be required to determine whether this 
co-occurrence is a happenstance of our subject selection or 
whether it points to a specific requirement of marked epigenetic 
reprogramming in the neoplastic transformation of thymic cells. 
The extent of DNMT3A and ASXLl mutations in thymoma should 
be explored, and treatments targeting epigenetic reprogramming 
should be considered in the future. 

MATERIALS AND METHODS 

Study consent 

The study participant provided written informed consent under a protocol 
approved by the institutional review board of Scripps. 

Tumor and normal sample 

The formalin-fixed, paraffin-embedded block containing the thymoma 
tumor was sectioned 1 1 times, with each section placed on a separate 
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slide. One section, at 4-|im thick, was stained and used to assess tumor 
area and cellularity. The boundaries of the tumor were marl<ed, and the 
tumor cellularity was estimated to be >85%. The additional 10 sections 
were cut at 10-|im thicknesses, and the designated (the boundaries were 
traced onto each slide using the stained and marked slide) tumor areas 
were extracted and isolated using the RecoverAII FFPE nucleic acid 
isolation kit (Ambion, Austin, TX, USA). Freshly drawn blood was obtained 
for the analysis of normal DNA. Genomic DNA was extracted from whole 
blood using the QIAamp system (Qiagen, Valencia, CA, USA). Sequence 
alignment and variant calling were performed against the reference 
human genome (NCBI 37/hg 19). 

Whole-exome sequencing, variant calling and filtration 
Whole-exome sequencing was pursued utilizing Agilent SureSelect exome 
hybridization followed by barcoding and sequencing of paired 100-bp 
reads on the equivalent of ~ 1 lane of an lllumina HiSeq2500 (lllumina, San 
Diego, CA, USA) instrument (samples were barcoded and distributed across 
multiple lanes). Read mapping and initial variant calling and quality 
filtration was performed using the standard BWA-GATK^^'^^ best practices 
variant quality score recalibration approach. A mean coverage of 157.3X 
and 98.7X was achieved for the tumor and normal sample with 96.1% and 
98.7% of the target exome covered by >10 reads, respectively. 

Somatic point mutations were detected using MUTECT.^^ Somatic indels 
were detected using SomaticlndelDetector.^^ Somatic variants were then 
annotated using the SG-ADVISER (genomics.scripps.edu/ADVISER/) and 
Cypher Genomics, Inc. (San Diego, CA, USA) systems for further filtration. 
Somatic variants present at >0.1% frequency in the HapMap,^^ 1000 
Genomes,^^ NHLBI GO Exome Sequencing Project,^^ or 385 genomes from 
the Scripps Wellderly population, as well as variants present in segmental 
duplication regions that are prone to produce false-positive variant calls 
were removed. Somatic indels were stringently filtered to remove indels 
within simple repeats, repeats identical to the inserted/deleted sequence 
or adjacent to homopolymer repeats of length 5 + . Nonsynonymous, 
in-frame, frameshift, nonsense or canonical splice-site donor/acceptor site 
variants were retained for further analysis. 

Low-pass whole-genome sequencing, structural variant calling 
and filtration 

Low-pass whole-genome sequencing was performed by barcoding and 
sequencing of paired 100-bp reads on the equivalent of ~1 lane of an 
lllumina HiSeq2500. A mean coverage of 7.4X and 7.7X was achieved for 
the tumor and normal sample with 92.0% and 95.5% of the genome 
covered by >3 reads, respectively. Copy number and structural variants 
were identified by SVDetect.^^ For copy number variants, read density was 
determined for 100-kb windows at 10-kb intervals. The significance of read 
density differences between the tumor and normal sample were 
determined using a binomial test with the normal sample determining 
the probability of read mapping per window in order to adjust for 
mapability and other window-specific effects. Significant windows (P-value 
threshold = 1 E-8) were stringently filtered to retain called copy number 
variations only where the ratio of tumor to normal reads, adjusted for the 
difference in total number of reads, was at least 1.5X greater for copy 
number gains or 0.5X lower for copy number losses. For structural variants, 
an initial set of raw candidate structural variants was called using SVDetect 
with at least seven reads supporting structural variant calls in the tumor 
sample and three reads supporting calls in the normal sample. Tumor 
sample candidate structural variants were then filtered to remove 
structural variants with support in the normal sample or mapping to 
centromeric regions. 

Circos plots were generated via RCircos.^^ Mutation plots were 
generated via adaption of cBioPortal visualization plots.^^ 
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